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Abstract 

The propagation of short life time particles with consequently broad mass width are discussed 
in the context of transport descriptions. In the first part some known properties of finite life time 
particles such as resonances are reviewed and discussed at the example of the />meson. Grave 
deficiencies in some of the transport treatment of broad resonances are disclosed and quantified. 
The second part addresses the derivation of transport equations which permit to account for the 
damping width of the particles. Baym's ^-derivable method is used to derive a self-consistent and 
conserving scheme, which fulfils detailed balance relations even in the case of particles with broad 
■ mass distributions. For this scheme a conserved energy-momentum tensor can be constructed. 

t— I | Furthermore, a kinetic entropy can be derived which besides the standard quasi-particle part also 

J> ' includes contributions from fluctuations. 

o : 
g ■ 

; 1 Introduction and Prospects 

^ l With the aim to describe the collision of two nuclei at intermediate or even high energies one is 
i ■ confronted with the fact that the dynamics has to include particles like the A 33 or p-meson resonances 
^ ■ with life-times of less than 2 fm/c or equivalently with damping rates above 100 MeV. Also the collision 
Q\ , rates deduced from presently used transport codes are comparable in magnitude, whereas typical mean 
Q ! kinetic energies as given by the temperature range between 70 to 150 MeV depending on beam energy. 

i \ Thus, the damping width of most of the constituents in the system can no longer be treated as a 
0*3 ■ perturbation. 

As a consequence the mass spectra of the particles in dense matter are no longer sharp delta 



functions but rather acquire a width due to collisions and decays. The corresponding quantum prop- 
agators G (Green's functions) are no longer the ones as in standard text books for fixed mass, but 
£3 | have to be folded over a spectral function A(e,p) of finite width. One thus comes to a picture which 
^ ■ unifies resonances which have already a decay width in vacuum with the "states" of particles in dense 
^ ■ matter, which obtain a width due to collisions (collisional broadening). The theoretical concepts for a 
^ . proper many body description in terms of a real time non equilibrium field theory have already been 
^ ! devised by Schwinger, Kadanoff, Baym and Keldysh |IJ in the early sixties. First investigations of the 
' quantum effects on the Boltzmann collision term were given Danielewicz [[|, the principal conceptual 
problems on the level of quantum field theory were investigated by Landsmann || , while applications 
which seriously include the finite width of the particles in transport descriptions were carried out only 
in recent times, e.g.p,|4|-|T0f For resonances, e.g. the A 33 -resonance, it was natural to consider broad 



mass distributions and ad hoc recipes have been invented to include this in transport simulation mod- 
els. However, many of these recipes are not correct as they violate some basic principles like detailed 
balance 0], and the description of resonances in dense matter has to be improved. 

In this talk the transport dynamics of short life time particles are reviewed and discussed. In the 
first part some known properties of resonances are presented. These concern the equilibrium and low 
density (virial) limits. Some example discussions are given for the di-lepton spectrum resulting from 
the decay of p- mesons in a dense nuclear environment, both in thermal equilibrium and in a quasi- free 
scattering process. On the basis of this some deficiencies of presently used transport codes for the 
treatment of broad resonances are disclosed and quantified. They affect the di-lepton spectra already 
on a qualitative level and signal that the low mass side is grossly underestimated in the respective 
calculations. This motivates the question discussed in the second part, namely, how to come to a 
self-consistent, conserving and thermodynamically consistent transport description of particles with 
finite mass width. The conceptual starting point will be a formulation within the real-time non- 
equilibrium field theory. The derivation is based on and generalizes Baym's ^-functional method 



TTJ. The first-order gradient approximation provides a set of coupled equations of time- irreversible 
generalized kinetic equations for the slowly varying space-time part of the phase-space distributions 
supplemented by retarded equations. The latter account for the fast micro-scale dynamics represented 
by the four-momentum part of the distributions. Functional methods permit to derive a conserved 
energy-momentum tensor which also includes corrections arising from fluctuations besides the standard 
quasi-particle terms. Memory effects |Tl2j{L4jl appearing in collision term diagrams of higher order as 
well as the formulation of a non-equilibrium kinetic entropy flow can also be addressed [[U . 



2 Preliminaries 

The standard text-book transition rate in terms of Fermi's golden rule, e.g. for the photon radiation 
from some initial state \i) with occupation to final states |/) 



W = E lf rn(l-n f ) 
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with occupation n w for the photon, is limited to the concept of asymptotic states. It is therefore 
inappropriate for problems which deal with particles of finite life time. One rather has to go to the 
"closed" diagram picture, where the same rate emerges as 

W= "{^"^)-" + (2) 

with now two types of vertices — and + for the time-ordered and the anti-time ordered parts of the 
square of the amplitude. Together with the orientation of the - — > and — -> propagator lines one obtains 
unique diagrammatic rules for the calculation of rates rather than amplitudes. The just mentioned 
propagator lines define the densities of occupied states or those of available states, respectively. There- 
fore all standard diagrammatic rules can be used again. One simply has to extend those rules to the 
two types of vertices with marks — and + and the corresponding 4 propagators, the usual time-ordered 
propagator — > between two — vertices, the anti-time-ordered one — > between two + vertices and 
the mixed — — ► or — i ones with fixed operator ordering (Wightman-functions) as densities of occupied 



and available states. For details I refer to the textbook of Lifshitz and Pitaevski ||15 
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Fig. 1: Closed real-time contour with two external points x,y on the contour. 

Equivalently the non-equilibrium theory can entirely be formulated on one special time contour, 
the so called closed time path Jl|, fig. 1, with the time argument running from some initial time to to 
infinity and back with external points placed on this contour, e.g., for the four different components 
of Green's functions or self energies. The special — h or H — components of the self energies define the 
gain and loss terms in transport problems, of. eq. (|2|) and eqs. below. 

The advantage of the formulation in terms of "correlation" diagrams, which no longer refer to 
amplitudes but directly to physical observables, like rates, is that now one is no longer restricted to 
the concept of asymptotic states. Rather all internal lines, also the ones which originally referred to 
the "in" or "out" states are now treated on equal footing. Therefore now one can deal with "states" 
which have a broad mass spectrum. The corresponding Wigner densities i— > or — i are then no longer 



on-shell ^-functions in energy (on-mass shell) but rather acquire a width, as we shall discuss in more 
detail. 

For slightly inhomogeneous and slowly evolving systems, the degrees of freedom can be subdivided 
into rapid and slow ones. Any kinetic approximation is essentially based on this assumption. Then for 
any two-point function F(x,y), one separates the variable £ = (£1 — t2,f[ — r^), which relates to the 
rapid and short-ranged microscopic processes, and the variable X = + fa, f[ + r^), which refers to 
slow and long-ranged collective motions. The Wigner transformation, i.e. the Fourier transformation 
in four-space difference £ = x — y to four-momentum p of the contour decomposed components of any 
two-point contour function 

F^(X-p) = J + where i,j e {-+} (3) 

leads to a (co- variant) four phase-space formulation of two-point functions. The Wigner transformation 
of Dyson's equation fll9|) in { — h} notation is straight forward. For details and the extensions to include 
the coupling to classical field equations we refer to ref. |16|| . 

Standard transport descriptions usually involve two approximation steps: (i) the gradient expan- 
sion for the slow degrees of freedom, as well as (ii) the quasi-particle approximation for rapid ones. 
We intend to avoid the latter approximation and will solely deal with the gradient approximation for 
slow collective motions by performing the gradient expansion of the coupled Dyson equations. This 
step indeed preserves all the invariances of the $ functional in a ^-derivable approximation. 

It is helpful to avoid all the imaginary factors inherent in the standard Green's function formula- 
tion and change to quantities which are real and positive either in the homogeneous or in appropriate 
coarse graining limits. They then have a straight physical interpretation analogously to the Boltzmann 
equation. We define 



F{X,p) = A(X,p)f(X,p) = T iG-+{X,p), 
F(X,p) = A(X,p)[l t f(X,p)\ = iG + -(X,p), 



with A(X, p) = -21m G R (X, p) = F±F (4) 



for the generalized Wigner functions F and F with the corresponding four phase space distribution 
functions f(X,p), the Fermi/Bose factors [1 =F f(X,p)] and spectral function A(X,p). According to 
the retarded relations between Green's functions G 13 , only two of these real functions are required for a 
complete dynamical description. Here and below upper signs relate to fermion quantities, whereas lower 



signs refer to boson quantities. As shown in ref. [16] mean fields and condensates, i.e. non- vanishing 



expectation values of one-point functions can also be included. 



3 Thermodynamic Equilibrium 

The thermodynamic equilibrium leads to a lot of simplifying relations among the kinetic quantities. 
All quantities become space-time independent. The Kubo-Martin-Schwinger condition determines the 
distribution functions to be of Fermi-Dirac or Bose-Einstein type, respectively 

f eq (X,p) = 1/ (exp ((p - n)/T) ± 1) . (5) 

Here \x is the chemical potential. The spectral function attains a form 

A (Y \ - r(p) -ih I r(p) = -2ImS R (p), 

eql ,P) M 2 (p) + r 2 (p)/4 Wlt \M(p) = M (p)-ReS R (p), M (p) = p% - pg(p). W 

This form is exact through the four- momentum p = (po,p) dependence of the retarded self-energy 
£ R (j>). Thereby Mo(p) = Pq — Pq(p) = is the free dispersion relation with k = 1 or 2 for the 
non-relativistic Schrodinger or the relativistic Klein-Gordon case, respectively. In the non-equilibrium 
case all quantities become functions of the space-time coordinates X and, of course, the distribution 
functions f(X,p) generally also depend on three momentum p. 



4 The Virial Limit 



Another simplifying case is provided by the low density limit, i.e. the virial limit. Since Beth-Uhlenbeck 
(1937) |17| it is known that the corrections to the level density are given by the asymptotic properties 
of binary scattering processes, i.e. in a partial wave decomposition by means of phase-shifts, see also 
[p"8| - p0|l . The reasoning can be summarized as follows. While for any pair the cm. motion remains 
unaltered the relative motion is affected by the mutual two-body interaction. Considering a large 
quantization volume of radius R and a partial wave of angular momentum j, the levels follow the 
quantization condition 




kR + 8j(E) 



(7) 



where Sj(E) is the phase-shift at relative energy E and n is an integer counting the levels. The kR 
term accounts for the free motion part. The corresponding corrections to both, the level density and 
thermodynamic partition sum, are given by 
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Fig. 2: 7T + 7r p-wave phase-shifts 
and scattering diagram. 



dn 
dE 



dn 



free 



dE 



+ 



2j + l d5j 



7T 



-Ei/T 



dE 
dn 



J dE 



(9) 



Since Z determines the equation of state, its low density limit 
is uniquely given by the scattering phase-shifts. The energy 
derivatives of the phase-shifts are also responsible for the time- 
delays discussed in ref. pT] and also for the virial corrections 
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to the Boltzmann collision term recently discussed in ref. 
The latter is directly connected to the B-term in our gener- 
alized kinetic equation fl2~0|). The advance of a phase-shift by 
a value range of 7r across a certain energy window adds one 
state to the level density and points towards an s-channel res- 
onance. An example is the p-meson in the p-wave 7r + 7r~ scat- 
tering channel, fig. 2. In cases, where the resonance couples 
to one asymptotic channel only, the corresponding phase-shifts 
relate to the vacuum spectral function Aj(p) of that resonance 
via* 



4|T ir 



r in (£)r out (£) 



(E* - E« R (E)Y 



F 2 

1 tot 



(E)/4 



(10) 



4 sin 2 S^E) = Aj(E,p = 0) T tot (E). (11) 



Here T inout is the corresponding T-matrix element. While relation ( |10D is correct also in the case where 
many channels couple to the same resonance, relation (O) only holds for the single channel case, 
where r in = r out = r tot . Relation (|TT|) illustrates that the vacuum spectral functions of resonances 
can almost model-independently be deduced from phase-shift information. In the case of the p-meson 
additional information is provided by the pion form factor. Also in the case of two channels coupling 
to a resonance the energy dependence of phase-shifts of the two scattering channels together with the 
inelasticity coefficient provide stringent constraints for the spectral function of the resonance 



" E is the relative cm. energy and correspondingly the momentum in A vanishes; 



1 for non-rel. particles; n = 2 



for relativistic bosons, where T{E)/2E equals the energy dependent decay width and E^E) 



V 2 



Re S(p). 



5 The p-meson in dense matter 



An an example I like to discuss the properties of the p-meson and the consequences for the decay into 
di-leptons. The exact production rate of di-leptons is given by the following formula 



dri 



e+e- 



dtdm 




f p (m,p,x,t) AJm,p,x,t) r" e+e "(m). (12) 



Here r pe+e (m) oc 1/m 2 is the mass- dependent electromagnetic decay rate of the p-meson into the 
di-electron channel. The phase-space distribution f p (m,p, x, t) and the spectral function A p (m,p,x,t) 
define the properties of the p-meson at space-time point x, t. Both quantities are in principle to be 
determined dynamically by an appropriate transport model. However till to-date the spectral functions 
are not treated dynamically in most of the present transport models. Rather one employs on-shell 
5-functions for all stable particles and spectral functions fixed to the vacuum shape for resonances. 

As an illustration the model case is discussed, where the p-meson just strongly couples to two 
channels: naturally the 7r + 7r~ channel and to the ttN <-* pN channels relevant at finite nuclear den- 
sities. The latter component is representative for all channels contributing to the so-called direct p 
in transport codes. For a first orientation the equilibrium properties are discussed. Admittedly by 
far more sophisticated and in parts unitary consistent equilibrium calculations have already be pre- 
sented in the literature, e.g. pi| - P^ ]. It is not the point to compete with them at this place. Rather 
we try to give a detailed analysis in simple terms with the aim to discuss the consequences for the 
implementation of such resonance processes into dynamical transport simulation codes. 

Both considered processes add to the total width of the p-meson 

T tot (m,p) = T p ^ 7T+7T -(m,p) +T p ^ NN -i(m,p), (13) 
and the equilibrium spectral function then results from the cuts of the two diagrams 

A p (m.p) = €3=( / ■ < =( 7 )=P = - Tp * + "~ + %^ . (11) 





( m 2_ m 2_ ReS ) +r t 2 ot /4 



In principle both diagrams have to be calculated by fully self consistent propagators, i.e. with cor- 
responding widths for all particles involved. This formidable task has not been done yet. Using 
micro-reversibility and the properties of thermal distributions the two terms in (]TJ]) contributing 
to the di-lepton yield fll2|) can indeed approximately be reformulated as the thermal average of a 
7r + 7r - — > p — > e + e~-annihilation process and a ttN — > pN — > e + e~A-scattering process, i.e. 

^ (U+U- Vmr cr(n + -K~ -> p -> e + e") + f w f N v wN a(nN -> pA^ -> e + e"A)^ (15) 

However, the important fact to be noticed is that in order to preserve unitarity the corresponding 
cross sections are no longer the free ones, as given by the vacuum decay width in the denominator, 
but rather involve the medium dependent total width (|13|). This illustrates in simple terms that rates 
of broad resonances can no longer simply be added in a perturbative way. Since it concerns a coupled 
channel problem there is a cross talk between the different channels to the extent that the common 
resonance propagator attains the total width arising from all partial widths feeding and depopulating 
the resonance. While a perturbative treatment with free cross sections in ( |T5D would enhance the yield 
at resonance, m = m p , if a channel is added, c.f. fig. 2 left part, the correct treatment (0) even inverts 
the trend and indeed depletes the yield at resonance, right part in fig. 2. Furthermore one sees that 
only the total yield involves the spectral function, while any partial cross section only refers to that 
partial term with the corresponding partial width in the numerator! Unfortunately so far all these 
facts have been ignored or even overlooked in the present transport treatment of broad resonances. 



Di-lepton rates from thermal p-mesons (T — 110 MeV) 
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Quasi-free ttN collisions 
and spectral function 
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Fig. 3: e + e~ rates (arb. units) as a function of the invariant pair 
mass m at T = 110 MeV from 7r + 7r~ annihilation (dotted line) and 
direct p- meson contribution (dashed line), the full line gives the sum 
of both contributions. Left part: using the free cross section recipe, 
i.e. with Tt 
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The calculation are done with 
and F p ^ 7rNN -i(m p )/2m p = 70 MeV. 



right part for the correct partial rates 
,„(m p )/2m p = 150 MeV 



Fig. 4: Fermi motion averaged ttN — > 
pN — * e + e~N cross sections at pion 
beam momenta of 1 and 1.3 GeV/c 
(dashed and full curve) as a function 
of invariant pair mass m. The dotted 
line gives the spectral function used 
here and in fig. 3. 



Compared to the spectral function (dotted line in fig. 4) both thermal components in fig. 3 show a 
significant enhancement on the low mass side and a strong depletion at high masses due to the thermal 
weight / oc exp(— Po/T) in the rate (|12D . A similar effect is seen in genuine non-equilibrium processes 
like the di-lepton yield resulting from Fermi-motion averaged itN — > pN scattering, fig. 4. The latter is 
representive for the first chance collision in a irA reaction and shows a behavior significantly different 
from that obtained in refs. 1301 . For orientation the sub-threshold conditions for the two beam momenta 



are given by the vertical lines in fig. 4! 

Much of the physics can already be discussed by observing that the partial widths are essentially 
given by the type of coupling (s or p-wave, I = 0, 1) and the phase space available for the decay 
channel. For point-like couplings and two-body phase space or approximately in the case of one light 
particle among only heavy ones in the decay channel (e.g. for ttNN -1 ) one finds 
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(16) 



In the 7T7T case the corresponding strength is approximately given by the vacuum decay, while it depends 
on the nuclear density in the pN <-> ttN case. The simple phase-space behavior fll6|) suggests that far 



away from all thresholds (m 2 3> s t hr) the ratio (7r?T-annihilation) / (direct p) of the two components 



for m > 500 MeV 
however 



This 
no such 



should result to a fairly smooth and almost constant function of m, e.g 
kinematical constraint is nicely confirmed by some calculations, e.g. of ref. 
behavior is seen in the direct p-mesons so far computed in refs. |3(| |31|[J]. 

For completeness and as a stimulus for improvements the discussed defects in some of presently 
used transport treatments of broad resonances (vector mesons in particular) are listed below. The last 
column gives an estimate of multiplicative factors needed to restore the defect (uir is the resonance 
mass and T* (between 70 and 120 MeV) is a typical slope parameter for the corresponding beam 
energy). Many of the points are known and trivial, and can trivially be implemented. They are just 
of kinematical origin. However the associated defects are by no means minor. Rather they ignore 
essential features of the dynamics with consequences already on the qualitative level and affect the 
spectra by far more than any of the currently discussed im-medium effects, e.g. of the p-meson. 

tin refs. [|30| , pl| the direct p component appears almost like the spectral function itself, i.e. untouched from any 
phase-space constraints which come in through the distributions f p (X,p). The latter favour low masses and deplete 
the high mass components! In fact rather than being almost constant the ratios (irn) / (direct p) exhibit an exponential 
behavior exp(— m/T*) for m > 500 MeV with T* between 70 - 110 MeV depending on beam energy, pointing towards 
a major deficiency in the account of phase-space constraints for the direct p-meson component in these calculations. 



List of defects in some of the transport codes 

[a] The differential mass information contained in the distribution functions 
f(X,p) = f(X,m,p) of resonances is ignored and only the integrated total 
number is evaluated as a function of space-time (direct p in refs. [[?(], |3l||). 

[b] Except for the ir + ir~ — > e + e~ case most resonance production cross sections 
are parametrized such that they vanishes for y/s values below the nominal 
threshold, e.g. below ttin + ni p in the case irN — > pN. This violates detailed 
balance, since broad resonances can decay for m < itlr. 

[c] In partial cross sections leading to a resonance the randomly chosen mass is 
normally selected according to the spectral function. This is not correct since 
the corresponding partial width in the numerator of ( |14| ) has to be considered. 

[d] Different partial cross sections are simply added without adjusting the total 
width in the resonance propagator accordingly. This violates unitarity. 

[e] The Monte Carlo implementation of selecting the random mass m of the 
resonance (item [c]) is sometimes falsely implemented, namely ignoring the 
kinetic phase-space of genuine multi-particle final state configurations, e.g. in 
ttN — > pN . Applies also to the A-resonance , e.g. for NN — > NA. 

[f] For the electromagnetic decay of vector mesons some authors use a mass in- 
dependent decay rate, e.g. r„_> e + e - /m = const., rather than that resulting 
from vector dominance and QED with r p ^ e + e - oc 1/m 2 . 



restoring factor 

exp(— (m — mji)/T*) 
(factor 10 or more at 
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6 ^-derivable approximations 



The preceding section has shown that one needs a transport scheme adapted for broad resonances. 
Besides the conservation laws it should comply with requirements of unitarity and detailed balance. 
A practical suggestion has been given in ref. [f|] in terms of cross section prescriptions. However this 
picture is tied to the concept of asymptotic states and therefore not well suited for the general case, 
in particular if more than one channel feeds into a broad resonance. Therefore we suggest to revive 
the so-called ^-derivable scheme, originally proposed by Baym [|ll[] on the basis of a formulation of 
the generating functional or partition sum given by Luttinger, Ward ||32|| , and later reformulated in 
terms of path-integrals P3| . This functional can be generalized to the real time case (for details see 
[ILSf ) with the diagrammatic representation^ 



ir{G} = ir {G } + ]r 



& GD 



(17) 
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Thereby the key quantity is the auxiliary functional $ given by two-particle irreducible vacuum dia- 
grams. It solely depends on fully re-summed, i.e. self consistently generated propagators G(x, y) (thick 
lines). The consistency is provided by the fact that $ is the generating functional for the re-summed 
self-energy E(x, y) via functional variation of $ with respect to any propagator G(y,x), i.e. 

- i£ = T 8i$/5iG. (18) 

The Dyson equations of motion directly follow from the stationarity condition of T (|T7D with respect 
to variations of G on the contour^ 



5T{G}/5G = 0, (Dyson eq.) 



(19) 

* counts the number of self-energy E-insertions in the ring diagrams, while for the closed diagram of $ the value 
ri\ counts the number of vertices building up the functional 

§an extension to include classical fields or condensates into the scheme is presented in ref. p6[ 



In graphical terms, the variation ([Tj|) with respect to G is realized by opening a propagator line in 
all diagrams of The resulting set of thus opened diagrams must then be that of proper skeleton 
diagrams of S in terms of full propagators, i.e. void of any self-energy insertions. As a consequence, 
the ^-diagrams have to be two-particle irreducible (label c2), i.e. they cannot be decomposed into two 
pieces by cutting two propagator lines. 

The clue is that truncating the auxiliary functional $ to a limited subset of diagrams leads to 
a self consistent, i.e closed, approximation scheme. Thereby the approximate forms of <|>( a PP r -) define 
effective theories, where $( a PP r -) serves as a generating functional for the approximate self-energies 
£( appr, )(a;, y) through relation (0), which then enter as driving terms for the Dyson equations of 
the different species in the system. As Baym [|ll[ has shown such a ^-derivable approximation is 



conserving for all conservation laws related to the global symmetries of the original theory and at the 
same time thermodynamically consistent. The latter automatically implies correct detailed balance 
relations between the various transport processes. For multicomponent systems it leads to a actio = 
reactio principle. This implies that the properties of one species are not changed by the interaction 
with other species without affecting the properties of the latter ones, too. The ^-derivable scheme 
offers a natural and consistent way to account for this principle. Some thermodynamic examples have 
been considered recently, e.g., for the interacting 7riVA system || and for a relativistic QED plasma 



7 Generalized Kinetic Equation 

In terms of the kinetic notation (H) and in the first gradient approximation the generalized kinetic 
equation for F takes the form 

VF(X,p) = B in (X,p) + C(X,p) (20) 
with the drift term determined from the "mass" function (c.f. (^)) 

M(X,p) = M (p) - Re Z R (X,p) (21) 

through the Poisson bracket T>F = {M, F}. The explicit form of the differential drift operator reads 

y dp" J x dX» dp^ d PfM \2p" rel. bosons. v ; 

The two other terms in (|20|) , B in (X, p) and C(X,p), are a fluctuation term and the collision term, 
respectively 

An = {r in , Re G R } , C(X,p) = r in (X,p)F(X,p) - T out (X,p)F(X,p), . (23) 
Here the reduced gain and loss rates and total width of the collision integral are 

T in (X,p) = T iX-+(X,p), T out (X,p) = iX+-(X,p), (24) 

r(x, P ) = -2im e r (x, p ) = r mt {x, P ) ± r ta (x, P ). (25) 



The combination opposite to ( |25D determines the fluctuations 

i(x, P ) = r in (x,p)Tr out (x, P ). (26) 

We need still one more equation, which in fact can be provided by the retarded Dyson equation. 
In first order gradient approximation the latter is completely solved algebraically 

nR = I ^ A(Xn) = r (^'P) (07) 

M{X,p)+iT{X,p)/2 1 ,P) M 2 (X,p) + r 2 (A,p)/4 1 } 



Canonical equal-time (anti) commutation relations for (fermionic) bosonic field operators provide the 
standard sum-rule for the spectral function. 



We now provide a physical interpretation of the various terms in the generalized kinetic equation 
fl20|). The drift term T>F on the l.h.s. of eq. fl20|) is the usual kinetic drift term including the 
corrections from the self-consistent field Re E fl into the convective transfer of real and also virtual 
particles. For the collision-less case C — B — 0, i.e. T>F = (Vlasov equation), the quasi-linear first 
order differential operator T> defines characteristic curves. They are the standard classical paths in the 
Vlasov case. Thereby the four-phase-space probability F(X,p) is conserved along these paths. The 
formulation in terms of a Poisson bracket in four dimensions implies a generalized Liouville theorem. 
For the collisional case both, the collision term C and the fluctuation term B change the phase-space 
probabilities of the "generalized" particles during their propagation along the "generalized" classical 
paths given by T>. We use the term "generalized" in order to emphasize that particles are no longer 
bound to their mass-shell, M = 0, during propagation due to the collision term, i.e. due decay, creation 
or scattering processes. 

The r.h.s. of eq. fl2"0|) specifies the collision term C in terms of gain and loss terms, which also 
can account for multi-particle processes. Since F includes a factor A, the C term further deviates from 
the standard Boltzmann-type form in as much that it is multiplied by the spectral function A, which 
accounts for the finite width of the particles. 

The additional Poisson-bracket term 



d _r F r> pr n\ ^-r 2 /4 Mr 

n in - y in , rve ur j _ ^ ^ u i in -h + ^ |i in , i out j- 



is special. It contains genuine contributions from the finite mass width of the particles and describes 
the response of the surrounding matter due to fluctuations. This can be seen from the conservation 



laws discussed below. In particular the first term in (28) gives rise to a back-flow component of 
the surrounding matter. It restores the Noether currents to be conserved rather than the intuitively 
expected sum of convective currents arising from the convective T>F terms in (|20|) . The second term 
of (|28| ) gives no contribution in the quasi-particle limit of small damping width limit and represents a 
specific off mass-shell response due to fluctuations, c.f. [^, [TJ|]. In the low density and quasi-particle 



limit the B m term provides the virial corrections to the Boltzmann collision term |22| . 



8 Conservations of the Current and Energy— Momentum 



The global symmetries of $ provide conservation laws such as the conservation of charge and energy- 
momentum. The corresponding Noether charge current and Noether energy-momentum tensor result 
to the following expressions, c.f. [TO, 



f{X) = 6 -Tr J ^Lv» (F(X,p) T F(X,p) 
&»(x) = ^Tr J 701^ (F(X,p) T F(X,pj) + gT (^PO " £ P °W 



(29) 



Here 



S mt (X) = (-C mt (X) 



<5$ 



5X(x) 



A=l 



£P Qt = -Tr 



/ 4^ [ Re sR ( ftf ) + Re ° R ( Fin T To ^ 



(27 



are the densities of the interaction energy and the potential energy, respectively. The first term of 
gpot com plies with quasi-particle expectations, namely mean potential times density, the second term 
displays the role of fluctuations / = r in =p r out in the potential energy density. This fluctuation term 
precisely arises form the 5-term in the kinetic eq. (|20D , discussed around eq. (^8|). It restores that the 
Noether expressions ( |2D| ) are indeed the exactly conserved quantities. In this compensation we see the 
essential role of the fluctuation term in the generalized kinetic equation. Dropping or approximating 
this term would spoil the conservation laws. Indeed, both expressions in (^) comply exactly with the 
generalized kinetic equation (EH), i.e. they are exact integrals of the generalized kinetic equations of 



motion within the ^-derivable scheme. Memory effects and the formulation of a kinetic entropy can 
likewise be addressed flU . 
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